Global trends in the interplane penetration depth of layered superconductors 
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We report on generic trends in the behavior of the interlayer penetration depth A c of several different 
classes of quasi two-dimensional superconductors including high-T c cuprates, S^RuCU, transition 
metal dichalcogenides and organic materials of the (BEDT — TTF)2X - series. Analysis of these 
trends reveals two distinct patterns in the scaling between the values of A c and the magnitude of 
the c-axis DC conductivity o~dc- one realized in the systems with a ground state formed out of 
well defined quasiparticles and the other seen in systems in which the quasiparticles are not well 
defined. The latter pattern is found primarily in under-doped cuprates and indicates a dramatic 
enhancement (factor ~ 10 2 ) of the energy scale Q.c associated with the formation of the condensate 
compared to the data for conventional materials. We discuss the implication of these results on the 
understanding of superconductivity in high-T c cuprates. 
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I. INTRODUCTION 

The formation of the superconducting condensate in el- 
emental metals and their alloys is well understood within 
the theory of Bardeen, Cooper and Schrieffer (BCS) in 
terms of a pairing instability in the ensemble of Fermi liq- 
uid (FL) quasiparticles. Applicability of the FL descrip- 
tion to high-T c cuprate superconductors is challenged by 
remarkable anomalies found in both the spin- and charge 
response of these compounds in the normal statefil. Be- 
cause quasiparticles are not well defined at T > T c in 
most cuprates it is natural to inquire into the distin- 
guishing characteristics of the superconducting conden- 
sate which appears to be built from entirely different 
" raw material" . Infrared spectroscopy is perfectly suited 
for this task, fndeed, the analysis of the optical con- 
stants in the far infrared (IR) unfolds the process of the 
formation of tha condensate <5(0)-peak in the dynami- 
cal conductivity^ and also gives insight into the single- 
particle excitations in the system both above and below 
T c . 

In this paper we focus on the inter-plane properties of 
high-T c superconductors. We will show that the distinc- 
tions in the behavior of the condensate in conventional 
superconductors and high-T c cuprates are most radical in 
the case of the c-axis inter-plane response. The analysis 
of the generic trends seen in the behavior of the c-axis 
condensate (correlation between the penetration depth 
A c and the DC conductivity ctdc) allows us to infer the 
energy scale ile associated with the development of the 
superfluid in the cuprates. This energy scale may dra- 
matically exceed the energy gap in systems lacking well 
defined quasiparticles at T > T c (primarily in under- 



doped cupratesQ). We discuss a connection between the 
magnitude of flc and the nature of the normal state re- 
sponse. 



II. EXPERIMENTAL PROCEDURES 

The response of the superconducting condensate can 
be investigated through the IR experiments probing the 
complex conductivity a(to) — <j\(oj) + ia^uS) of a super- 
conductor. At T < T c the real part of the conductivity 
can be written as: 
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The delta-peak term represents the response of the con- 
densate with the superfluid density p s = Aim. s e /m* pro- 
portional to the concentration of superconducting carri- 
ers n s and inversely proportional to their effective mass 
to* . The second term on the right-hand side of the Eq. |l| 
is usually referred to as the regular component and repre- 
sents the conductivity that is NOT due to the supercon- 
ducting carriers. It may include conductivity due to un- 
paired carriers at T < T c at finite frequencies, phonons, 
interband transitions, magnons, etc. Commonly, the con- 
densate stiffness is characterized through the penetration 
depth A = c/y/pj, the notation we will use in this paper. 

In order to discuss several techniques that can be ex- 
ploited to determine the interlayer penetration depth 
of an anisotropic superconductor we turn to our data 
for La 1 . 83 Sr .i 7 CuO 4 (La214) with T c ~ 36 K (Fig. @). 
Large single crystals were grown using traveling-solvent 
floating-zone techniqueQ and were carefully annealed to 
remove excess oxygen. The crystallographic axes were 
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determined by Laue diffraction and the samples were 
then cut into platelets with the ac planes parallel to 
the wide face. The error in the axes directions is less 
than 1°. Near-normal- incidence reflectance measure- 
ments were performed at UCSD in the frequency range 
between 10-48,000 cm" 1 (1 meV - 6 eV). The com- 
plex conductivity o~(lu) and complex dielectric function 
e(u>) — e±(ou) + 1^2(^2) were inferred from R(lu) using 
Kramers-Kronig (KK) analysis. The low and high fre- 
quency extrapolations have negligible effect on the data 
in the measured frequency interval. 

Below we outline common analysis techniques used to 
determine the penetration depth from the results of IR 
studies. 

1) Raw c-axis reflectance of most high-T c supercon- 
ductors at T<CT C exhibits a sharp plasma edge. In the 
case of Lai.83Sr .i7CuO4 this feature is located at ~85 
cm -1 (Fig. [|, panel A). This behavior is in contrast to 
the featureless normal state reflectance. The position of 
the plasma edge is determined by the screened plasma 
frequency tu p , from which the penetration depth can be 
obtained as c 2 /A 2 = oJ p 2 e 00 (Ref. ||). In the latter equa- 
tion is the real part of the dielectric constant €i(w) at 
frequencies above the plasma edge. The numerical value 
of £oo is somewhat ambiguous and introduces error in 
the result for A c . This technique has been employed in 
Ref. §-|. 

2) In a BCS superconductor the formation of the con- 
densate is adequately described with the Ferrel-Glover- 
Tinkham (FGT) sum rule: 

Ps = t= r°KH-<rr>)]^ (2) 
A Jo+ 

where erf (u) is the normal state conductivity at T c and 
the upper integration limit Qc is of the order of the 
gap energy. The upper cut-off issue for cuprates will 
be discussed in detail below. According to this sum 
rule the area "missing" from the normal state conduc- 
tivity (shaded region in Fig. [I] B) is recovered under 
the <5(0)-peak. This technique may somewhat underes- 
timate the magnitude of A c because, at least in under- 
doped cuprates the superfluid density is accumulated 
from a broad_ea£rgy region significantly exceeding the 
gap energyoi§E3Eil. This method has been used for the 
analysis of the penetration depth in Ref. [l2| [l4|. 

3) Finally, the most commonly used method of extract- 
ing A c is based on the examination of the imaginary part 
of the complex optical conductivity. By KK transforma- 
tion, the 5-peak at lu = in the real part of the optical 
conductivity implies that the imaginary part has the form 
<T2(w) = c 2 /(47rwA 2 ). Therefore the magnitude of A c can 
be estimated from lu x 02(0;) in the limit of lu — > 0. (the 
gray line in panel C of Fig. [j] or the dotted line in panel 
D) (Ref. 

While method 3 is very well suited to quantify the 
magnitude of the penetration depth, this technique also 
may introduce systematic errors. Strictly speaking the 



relation 02(0;) = c 2 /(47rwA 2 ) is valid only if er[ e£, (w) = 0. 
Typically, this is not the case in high-T c superconductors 
which all show residual absorption in the far-IR conduc- 
tivity. This absorption may be (in part) connected with, 
d-wave symmetry of the order parameter in cupratesu 
leading to gapless behavior at any finite temperature. 
Data displayed in Fig. |l| B clearly shows non-vanishing 
IR conductivity down to the lowest T and lu. A finite 
regular contribution to a\ (lu) implies a finite contribu- 
tion to <7 2 ( w )- Owing to this contribution the spectra of 
02(0;) acquire a complicated frequency dependence that 
may significantly differ from the l/u> form (Fig. [l]C, D). 
Moreover, the magnitude of the penetration depth ex- 
tracted from such a spectrum is likely to be underesti- 
mated, even if the product 02 (to) x lu is taken at the 
lowest experimentally accessible frequencies. 

Systematic errors in the magnitude of A connected with 
o\ e9 (to) > can be eliminated using the following pro- 
cedure. The intrinsic value of the penetration depth can 
still be determined from 02 (to), if the imaginary part of 
the conductivity is corrected by a r 2 e9 (tu) characterizing 
all screening effects that are not due to superconducting 
carriers at T < T c : 

c 2 

ff3 ( w )-a 2 r » = i _ x -. (3) 

To determine a 2 (to) we employ a KK-like transforma- 
tion: 

reg, s 2w f°° q[ £g (uf) , 

C 2 M = / — 2 d0J ■ 4 

7T J 0+ LU U - LU 

The result of the application of the correction procedure 
for the imaginary part of the conductivity is presented 
in Fig. [I] D. It appears that after subtraction of <j^ 9 (lu) 
the remaining contribution to the conductivity reveals 
al/u behavior over an extended frequency region, sup- 
porting the soundness of the procedure proposed here. 
We emphasize again that no other correction procedure 
besides that described by Eqs. || and |^ has been used. 
In the case of Lai.g3Sro.i7Cu04 the latter procedure has 
lead only to minor correction of the absolute value of A c 
(~ 18%). That is because the absolute value of er[ es (w) is 
relatively small and is constant throughout far-IR (Fig. |l| 
B). However, such a correction can be much more signifi- 
cant for the ovcrdoped samples which often show stronger 
Drude-like contribution in o\ e9 (lu) spectra. Panel D also 
shows a frequently used approximation to the method we 
have just outlined: instead of subtracting a r 2 eg (to), one 
subtracts <j2(to,T c ) from 02(00, T -C T c ). The resulting 
curve looks somewhat better than the uncorrected one, 
but still yields an enhanced value of lu x 02 (w) in the limit 
of lu -> 0. 



III. UNIVERSAL C-AXIS PLOT 

The c-axis penetration depth in a layered swpeFepn- 
ductor can be determined from IR experiments0Enl3 as 
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described in the previous section. In addition, several 
other experimental techniques including -magnetization 
measurement — c3, microwave absorptior£3~E3, and vor- 
tex imagingoEi can be used to determine the magnitude 
of A c . Regardless of the method employed, the inter- 
layer penetration depth in several families of cuprates 
reveals a universal scaling-behavior with the magnitude 
of <r DC {T = T c ) (Fig. the absolute value of A c is 

systematically suppressed with the increase of the nor- 
mal state conductivity^. The scaling is obeyed primar- 
ily in under-doped cuprates (blue symbols in Fig. ^J) . 
The deviations from the scaling are also systematic and 
are most prominent in over-doped phases (red symbols 
in Fig. ||). Such deviations are a direct consequence of a 
well-established fact: on the over-doped side of the phase 
diagram a^c increases whe.re.fi S.JL is either unchanged or 



may show a minor increase 

We find a similar scaling pattern between A c and o\dc 
in other classes of layered superconductors, including 
organic materials, transition metal dichalcogenides and 
Sr 2 Ru04 (Fig. ||). While the non-cuprate data set is 
not nearly as dense, the key trend is analogous to the 
one found for cuprates. The slope of the A c — croc de- 
pendence is also close for both cuprates and non-cuprate 
materials. The principal difference is that the cuprates 
universal line is shifted down by approximately one or- 
der of magnitude in A c . The latter result shows that the 
superfluid density (cx 1/A 2 ) is significantly enhanced in 
under-doped cuprates compared to non-cuprate materi- 
als with the same DC conductivity. 

Possible origins of the A c — a pa correlation were re- 
cently discussed in the literature^! A plausible qualita- 
tive account of this effect can be based on the FGT sum 
rule, Eq. H For a dirty limit superconductor (oj) w 
a DC: and Eq. ^| can be approximated as: 



2A<t DC - 



(5) 



Such an approximation is possible because within the 
BCS model the energy scale flc from which the conden- 
sate is collected is of the order of magnitude of the gap: 
f2c ~ 2A ~ 3 — 5kT c . A connection between 1/A 2 , (Jdc, 
and 2A is illustrated in the inset of Fig. |[ In the dirty 
limit the magnitude of odc sets the amount of spectral 
weight available in the normal state conductivity whereas 
the magnitude of flc — 2A defines the fraction of this 
weight which is transferred into condensate at T < T c . 
Therefore, the magnitude of A c can be expected to sys- 
tematically decrease with the enhancement of the DC 
conductivity, in accord with the FGT sum rule. Notably, 
an approximate form (Eq. o) yields the A c — croc scaling 
with the power law a = 1/2 which is close to a = 0.59 
seen in Fig. |[ 

The strong condensate density in the cuprates can be 
understood in terms of the dramatic enhancement of the 
energy scale flc over the magnitude of the energy gap. 
This can be seen through a comparison of the universal 



scaling patterns observed for cuprates and of a similar 
pattern detected for non-cuprate superconductors. The 
energy scale associated with the condensate formation 
for materials on the upper line, which for most conven- 
tional materials in Fig. |^ is close to estimates of the gap, 
is of the order of 1-3 meV. In Sr2Ru04 for example 2A 
= 2.2 meV based on Andreev reflection measurements.!^ 
If Eq. H in the form p s = c 2 /A 2 « ftc<JDC appropriate 
for cuprates, is employed to describe the difference be- 
tween the upper and the lower lines in Fig. |^, then one 
can conclude that the corresponding scale for undcrdoped 
cuprates is ~ 100 times greater, i.e. of the order 0.1 - 0.3 
eV. This assessment of £lc is supported by the explicit 
sum rule analysis for several cupratesBO and also makes 
£lc the largest energy scale in the problem of cuprate 
superconductivity^!! 

Data points in Fig. || for overdoped materials support 
the notion that the A c — <jdc plot provides means to learn 
about the energy scale associated with the condensate 
formation. Deflection of the over-doped cuprates from 
the universal line implies that ftc is gradually suppressed 
with the increased carrier density. This trend is common 
for Tl 2 Ba 2 Cu0 6+5 (T12201), La214 and YBa 2 Cu 3 7 _ 5 
(YBCO) materials (see Fig. |[). Integration of the con- 
ductivity for all these overdoped materials shows that the 
FGT sum rule is exhausted at energies as low as 0.08 eV 
(Ref. 00). 

In BCS superconductors flc is related to 2A, and 
therefore to T c . In cuprates we find no obvious connec- 
tion between the broad energy scale £lc and the critical 
temperature T c . While scaling of A c by the magnitude, af 
T c does reduce the "scattering" of the data pointsEjEj, 
the two distinct A c — <jdc patterns persist even if such 
scaling is implemented. Similarly, the difference between 
the two lines in Fig. ^ cannot be accounted for by dif- 
ferences in T c . In particular, the critical temperature of 
strongly under-doped La214 materials is nearly the same 
as that of the several ET-compounds (~ 12 — 15 K). Nev- 
ertheless, the penetration depth is dramatically enhanced 
in the latter systems. 



IV. IN-PLANE QUASIPARTICLES AND 
INTER-PLANE TRANSPORT 

A quick inspection of the materials in Fig. |2| suggests 
that the ~ 3 meV scale (top line) is observed in systems 
in which superconductivity emerges out of a normal state 
with well defined quasiparticles, whereas the enhanced 
value of Qc — 0.3 eV is found in undcrdoped cuprates for 
which the quasiparticle concept may not apply (bottom 
line). The experiments which in our opinion are most rel- 
evant to this classification include quantum oscillations 
of the low-T inter-layer resistivity (and of other quanti- 
ties) in high magnetic fieldscJ. Quantum oscillations can 
be viewed as a direct testimony of long-lived quasiparti- 
cles capable of propagating coherently between the lay- 
ers. Indeed, quantum oscillations have been observed in 
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2D organic superconductors0@, 2H-NbSe 2 (Ref. [l3|) and 
Sr 2 Ru04 (Ref. [f4|). On the contrary, quantum oscilla- 
tions have never been reported for under-doped cuprates. 
The lack of coherence in the c-axis transport in these ma- 
terials indicates that the ground state of cuprates may be 
fundamentally different. 

Signatures of coherent and incoherent behavior can 
also be recognized in the spectra of the c-axis conduc- 
tivity. A hallmark of a coherent response is the Drude 
peak seen in <7i(w) of metals. Notably, a similar feature 
has never been found in underdoped compounds (form- 
ing the lower line in Fig. ^). The electronic contribu- 
tion to cti(w) in these materials is usually structureless 
which is commonly associated with the incoherent (dif- 
fusive) motion of charge carriers across the planes. On 
the contrary, many materials that belong to the upper 
line in Fig. g demonstrate a familiar Drude-like behav- 
ior. This kind of behavior has been found in Sr2Ru04 
(Ref. ^BJ) and is also shown in our datao for inter-plane 
response of 2H-NbSe 2 (Fig. || top-right panel). In both 
cases, the width of the peak decreases at low tempera- 
tures, which is characteristic of the response of ordinary 
metalscll. As for the over-doped cuprates (located in a 
cross-over region between the two lines in Fig. ^|) their 
conductivity is indicative of the formation of the Drudc- 
like peak (see for example a c {uo) for YBa2Cri307; Fig. [|, 
top-middle panel), which is becoming more pronounced 
with increased carrier densityo. 

Analysis of the anisotropic carrier dynamics in several 
layered superconductors indicates that the degree of co- 
herence in the interplane transport may be related to 
the strength of inelastic scattering within the conducting 
planes. The bottom panels in Fig. ^ show the in — plane 
scattering rate (inverse lifetime) l/r a b(u;) (Ref. [49|) for 
the layered compounds corresponding to the top three 
panelse3. In all these systems l/T a ;,(a->) oc u> ewer an 
extended frequency interval (up to 3,000 cm _1 )EJ. An 
important feature of the data displayed in Fig. is that 
as doping is increased from underdoped YBa2Cu306.6 
to optimally doped YBa2Cu30g.95 the absolute values of 
1/Ta&(w) decrease. j-A similar trend is observed in other 
cuprate familiesE3~E3. The shaded regions in Fig. || repre- 
sent Landau Fermi liquid (LFL) regime, where the quasi- 
particles are well defined, i.e. the magnitude of the scat- 
tering rate is smaller than energy (1/t(cj) < lu). In 2H- 
NbSe 2 1/T a b(uj, 10K) is in the LFL regime over the entire 
frequency interval displayed in Fig. ||. However this is not 
the case for the two cuprates discussed. We believe that 
these differences in absolute values may have a profound 
effect on the interplane transport. In 2H-NbSe2 where 
the in-plane quasiparticles are well defined the interplane 
transport is also coherent, and is characterized by a nar- 
row Drude-like mode whose width decreases with tem- 
perature (Fig. ||, top-right panel). On the other hand, 
in YBa2Cu306.6> which is lacking well defined quasipar- 
ticles, the interplane transport is incoherent, with ui(uj) 
being dominated by optical phonons (Fig. [SI top-right). 
As for the over-doped YBa2Cu30y (Fig. ||[ bottom- 



middle) the optical conductivity of this compound is in 
between these two opposite limits. Figure |3j therefore 
supports the notion that long-lived in — plane quasipar- 
ticles may be one of the necessary prerequisites for co- 
herent out — of — plane transport. 

V. GLOBAL TRENDS IN LAYERED 
SUPERCONDUCTORS 

To summarize the experimental results reported in this 
work, we wish to stress the following points: i) two dis- 
tinct patterns in A c — ctdc correlation (Fig. ^) are in- 
dicative of a dramatic difference (~ 10 2 ) in the energy 
scale flc from which the interlayer condensate is col- 
lected; ii) the pattern with the typical energy scale of 
£lc — 3 meV is realized in the materials with the coher- 
ent transport between the planes, whereas the one with 
£lc — 300 meV is found in underdoped cuprate super- 
conductors with an incoherent response; Hi) over-doped 
cuprates reveal a cross-over between the two behaviors; 
iv) coherence in the interlayer transport correlates with 
the strength of inelastic scattering within the conducting 
planes (Fig. ||). 

These results allow us to draw several conclusions re- 
garding features of the superconducting condensate in 
different layered systems: 

• The symmetry of the order parameter seems to be un- 
related to trends seen in the c-axis condensate response. 
Indeed, the upper line in Fig. |^ is formed by s-wave tran- 
sition metal dichalcogenides, p-wave Sr2Ru04 and or- 
ganic materiaLs-ibr which both s- and d-wave states have 
been proposedEll, while d-wave high-T c materials form the 
lower line and the crossover region between the lines. 

• Electrodynamics of the systems on the top line at 
T <C T c is determined by the magnitude of the gap (and 
hence by T c ), in general agreement with the BCS theory. 
It is therefore hardly surprising that the trend initiated 
by 2D superconductors is also followed in 1-dimensional 
organic conductors, as well as by more conventional sys- 
tems such as Nb Josephson junctions, bulk Nb and Pb 
or amorphous aMoi- x Ge x (see Fig. ||). 

• While the pseudogap state has been shown to be re- 
sponsible for the anomalous supcrfluid response of the 
underdoped cupratesaO, the characteristic pseudogap 
temperature T* = 90 — 350 K is still much lower than 
our estimate of fic for these materials (0.1 - 0.3 eV, i.e. 
1,000 - 3,000 K). 

• Unlike BCS superconductors where T c is determined 
by 2 A and therefore by f2cs the critical temperature T c 
in cuprates correlates with neither 2 A nor Qc- 

In conclusion, analyzing a large amount of experimen- 
tal data, we found two distinctly different patterns in 
A c — o-£,c scaling in layered superconductors. Based on 
the universal c-axis plot, we inferred a broad energy scale 
fl c relevant for pair formation in underdoped cuprates. 
This result is consistent with the idea that the supercon- 
ducting transition in the cuprates is driven by a lowering 
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of the electronic kinetic energyEl. We argue that the ap- 
pearance of such an energy scale is fundamentally related 
to the incoherent c-axis transport, which on the other 
hand may be related to poorly defined in-plane quasi- 
particles. A quantitative account of the distinct energy 
scales associated with the condensate is a challenge for 
models attempting to solve the puzzle of cuprate super- 
conductivity. 
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FIG. 1. Interlayer response of LSCO single crystals with 
T c =36 K: reflectance R(w) (panel A); real and imaginary 
parts of the conductivity (panels B and C) and the product 
02 (w) x uj (panel D). The c-axis penetration depth can be de- 
termined from the IR data using several different techniques: 
from the position of the plasma minimum in R(w), from inte- 
grating the difference between the <j\(lu,T c ) and o~i(lu, 10K) 
(Eq. □); and from examining the frequency dependence of the 
oi{u), WK) x uj. Advantages and deficiencies of these methods 
are analyzed in Section II. The latter approach may under- 
estimate the magnitude of A c because of the screening ef- 
fects associated with the response of unpaired charge carriers 
at T<CT C . We employed a Kramers-Kronig transformation 
(Eqs. ^ and|ij) to correct for this effect (solid line in panel D). 



FIG. 2. The c-axis penetration depth A C (T = OK) as a 
function of the c-axis DC conductivity <jdc(T c ). We find 
two distinct patterns of A c — <jdc scaling. Cuprate super- 
conductors exhibit much shorter penetration depths than 
non-cuprates materials with the same (jdc{T c ). This re- 
sult implies a dramatic enhancement of the energy scale Qc 
from which the condensate is collected as described in the 
text. The superconducting transition temperature T c has 
not been found to be relevant to the A c — ctdc scaling. Data 
points: YBCO (Ref.|F|,0), overdoped YBCOJRef. jl^-fLs!) 
La214 n (^eL_|,^|Jl|,|§; 
(Ref. " 
(Ref. |8l 



HgBa 2 Cu 2 04 (Ref. ^J5J), T12201 
IEJ35|), Bi 2 Sr 2 CaCu20 8 (Ref. @||), Nd 2 - !C Ce : < : Cu04 
1 27). Blue points - underdoped (UD), green - opti- 
mally doped (OpD); red - overdoped (OD). Transition metal 
dichalcogCTides^Ref. Fi[j2l]0,|o|-^|), (ET) 2 X compounds 



(Ref. ^]2j,^,^M63j-|6j), (TMTSFUCIO4 (Ref. |6j,|S7|), 
Sr 2 Ru0 4 (Ref. g6|||), niobium (Ref. EjEil ), lead (Ref. ||), 
niobium Josephson junctions (Ref. |39D and aMoi_ x Gei 
(Ref. [70]). Inset: in a conventional dirty limit superconduc- 
tor the spectral weight of the superconducting condensate 
(given by 1/A 2 ) is collected primarily from the energy gap 
region. The total normal weight is preset by magnitude of 
<jdc whereas the product of 2A x odc quantifies the fraction 
of the weight that condenses. 



FIG. 3. Examples of the interplane transport for lay- 
ered superconductors. Top panels show the out-of-plane 
optical conductivity o c (uj), the bottom panels the corre- 
sponding in-plane scattering rate 1/T a b(u>). The observa- 
tion of the Drude-like feature in the interplane optical con- 
ductivity of the dichalcogenide 2H-NbSe2 (top-right panel) 
is consistent with magnetoresistance measurements that re- 
vealed evidence for well-behaved quasiparticles. In con- 
trast the conductivity of underdoped YBa2Cua06.6 material 
(top-left panel) gives no signs of coherent response. Over- 
doped cuprates show the emergence of a Drude-like fea- 
ture (top-middle panel) and also occupy an intermediate 
position between the two lines in Fig. ^. Experimental 
data: YBa 2 Cu 3 6 .6 (Ref. |§fn]), YBa 2 Cu 3 07 (Ref. |l|@, 
YBaaCuaOe.gs (Ref. B§) and 2H-NbSe 2 (Ref. S). 
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